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Abstract 

We investigate LIBOR-based derivatives using a parsimonious field theory in- 
terest rate model capable of instilling imperfect correlation between different ma- 
turities. Delta and Gamma hedge parameters are derived for LIBOR Caps against 
fluctuations in underlying forward rates. An empirical illustration of our method- 
ology is conducted to demonstrate the influence of correlation on the hedging of 
interest rate risk. 
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1 Introduction 



LIBOR-based derivatives such as Caps and Floors are important financial contracts in- 
volving a sequence of quarterly payments ranging from one to ten years. Consequently, 
pricing and hedging such derivatives requires the modeling of multiple LIBOR rates. 

In an economy where LIBOR rates are perfectly correlated across different maturities, 
a single volatility function is sufficient. However, non-parallel movements in the LIBOR 
term structure introduce an important complication. To reduce the number of necessary 
inputs, volatility parameters within certain time intervals are often assumed to be iden- 
tical. However, this assumption represents a serious compromise, and longer maturity 
options still require a large number of volatility parameters even after such aggregation. 

In light of this issue, we utilize field theory models introduced by Baaquie pQ to 
instill imperfect correlation between LIBOR maturities as a parsimonious alternative to 
the existing theory. We derive the corresponding hedge parameters for LIBOR Caplets 
for applications to risk management. We then demonstrate the ease with which our 
formulation is implemented and the implications of correlation on the hedge parameters. 

Hedge parameters that minimize the risk associated with a finite number of random 
fluctuations in the forward interest rates is provided in Baaquie, Srikant, and Warachka 
[2]. Previously, field theory research has focused on applications involving traditional 
Heath, Jarrow, and Morton ^Hj forward interest rates, and on the pricing of LIBOR- 
based derivatives as is Baaquie J3]. This paper extends the concept of stochastic Delta 
hedging developed in Baaquie [I] to the hedging of LIBOR derivatives. 

The remainder of this paper begins with a review of the field theory model for pricing 
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LIBOR derivatives. Section 3 then investigates their corresponding hedge parameters, 
while Section 4 details their empirical implementation. The conclusion follows in Section 
5. 



2 Field Theory Model 

The introduction of imperfect correlation between all underlying LIBOR rates is ac- 
complished by the specification of a propagator for interest rate dynamics. In terms of 
notation, L(t, T) denotes the LIBOR rate at the current time t between time T and T + £ 
in the future where I = 1/4 year denotes the standard 3- month time interval between 
payoffs. 

Since forward rates are the basis for LIBOR rates, we first detail the Lagrangian 
underlying the evolution of forward rates. Let A(t, x) be a two dimensional field driving 
the evolution of forward rates f(t,x) through time 

= a(t,x) + a(t,x)A(t,x) (1) 

where a(t,x) and a(t, x) denote their volatility and drift velocity respectively. 

Following Baaquie and Bouchaud [3], the Lagrangian of the field is defined by three 
parameters, namely /i A and r\. 

Definition 2.1 The Lagrangian which describes the evolution of instantaneous forward 
rates equals 



where psychological future time is defined by z = (x — t) n . 



(2) 



The Lagrangian in Definition 12.11 contains a squared Laplacian term that describes 
the stiffness of the forward rate curve. Baaquie and Bouchaud j3] demonstrate that 
this formulation is empirically able to account for the phenomenology of interest rate 
dynamics. Ultimately, our pricing formulae for Caps and Floors stems from a volatility 
function and correlation parameters fi, X and r\ contained in the propagator, as well as 
the initial term structure. 

These forward rate dynamics are ultimately invoked for the pricing of Caps and Floors 
after expressing derivatives on interest rates in terms of their counterparts on bonds. 



2.1 LIBOR Dynamics 

The following is the relationship between the forward interest rates and the LIBOR term 
structure 

„/y +e dxf(t,x) _ 1 

Ht,T) = . (3) 

In the original Heath, Jarrow, and Morton model ^U], the martingale measure is 
defined by discounting Treasury Bonds denoted B(t,T) by the money market account 
R(t, t*), defined as 

R(t,U) = ef* Ur ® dt , (4) 

for the spot rate of interest denoted r(i). In contrast, in this paper all computations are 
carried out using the LIBOR measure for which LIBOR rates evolve as martingales. In 
other words, for t* > t 

L(t,T n ) = E L [L(U,T n )] . (5) 
Following the material in Baaquie |3], the drift ai{t, x) that corresponds to the LIBOR 
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martingale condition is given by 

ai{t,x) = —a(t,x) / dx D(x,x ;t)a(t,x ) ; T n < x < T n+ £ . (6) 

As proved in Baaquie 4j, a money market numeraire entails more complex calculations 
but arrives at identical prices if one instead uses the LIBOR measure. For the remainder 
of this paper, the subscript of L is suppressed with all expectations performed under the 
LIBOR measure. 

2.2 Pricing an Individual Cap let 

The existing literature justifies the Black model for pricing Caps and Floors by modifying 
risk neutral Heath, Jarrow, and Morton ^0] forward rates to yield LIBOR dynamics under 
the forward measure. Brace, Gatarek, and Musiela [7j is the seminal paper in this area, 
with additional details found in Musiela and Rutkowski [T2] . 

We review the field theory pricing formula for a Caplet for both a general volatility 
function a(t, T) and propagator D(x, x ; t) underlying risk neutral forward rates. Denote 
the principal amount of the Cap as V. If the Caplet is exercised at time T, the payment 
is made in arrears at time T + I. Hence the payoff function at time T + I is given by 

g(T + £)=iV(L(T,T)-K) + (7) 

where K denotes the strike rate of the Caplet. Note that before discounting the payoff 
at time T, we first discount from T + i back to time T . The entire expression for the 
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Caplet price is given by 

Caplet(t,T) = B(t,T)E m [B(T,T + £)g(T + £)] 



(8) 



V_ 
X 



B(t,T)E m 



fjC - e~ dxf(T,x) 



(9) 



according to equation (j^J and for X = l/(l + £K). Observe that invoking the forward 
measure involves multiplying by the bond B(t,T) with the only random forward rate 
term structure from T to T + I. Then, 



/+oo 
dG^(G,T,T + £){X -e 
-oo 



(10) 



where, as the derivation in Baaquie |T], \J/(C, T,T + £) equals 

1 



V 
X 



B(t,T). 



1 



27rg 2 (T - t) 6XP | 2g 2 (T-t) 



a - I T+l dxf(t,x) - g2(T o ^ ) } .(ii) 



The above result leads to the next proposition for Caplet pricing. 



Proposition 2.1 The price of a Caplet with strike K which matures at time T equals 



Caplet{t,T,T + 



V_ 
X 



B(t,T) [XN(d+) - FN{dL)] 



for X = , B(t,T) = 

i+iL(t,t,T) ' an d following definitions 



F 

d± 
q 2 



1 



l+£L(t,T) 
1 



qVT^t 
1 



T-t 



/' pT+t 

dt I dxdx a(t, x)D(x, x ;t)a(t, x ) 



(12) 



(13) 



Observe that the propagator for forward rates are elements of the Caplet price. The 
price of an at-the-money Caplet is then defined for X = F, which yields d± = J ^ q ^ T ~ t 



2 ' 
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implying an associated price of 

Caplet(t,T,T + £) = VB{t,T)[N{-dJ) - N{-d+)} (14) 



VB(t,T) 



N 



(15) 



3 Hedging 

This section details the implications of our field theory model on hedging LIBOR deriva- 
tives. The impact of correlation is examined in the context of the residual variance and 
the Delta hedge parameter for a portfolio. In particular, the more practical Stochastic 
Delta hedging technique is given in Subsection 13.21 



A portfolio II(t) composed of a Cap^t^T) 1 and N LIBOR futures contracts, with 
the futures chosen to ensure fluctuations in the value of the portfolio are minimized, is 
studied. This portfolio equals 

JV 

II(t) = Cap(t, U,T) + J2 ni{t)T{t, %) , (16) 

i=l 

where rii(t) represents the hedge parameter for the i th futures contract included in the 
portfolio. The LIBOR futures and Cap prices are denoted by 

•F(t,T;) = V[l-£L(t,T$ (17) 
Ca p(M.,T) = ^,,r)/ + °°^£=e-i(^"'-'-»-*) a ,X-e-) + (18, 

where for the midcurve Cap, we have 

l r u ~ f tt+£ 

ql = / dt J dxdx a(t, x)D(x, x ;t)a(t, x ) . (19) 



t Jt, 



lr This is a more general expression for a Cap referred to as the midcurve Cap. 
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N 



From equation (JTBJ), we have 

II(t) = Cap(t, t*,T) + VJ2 rn(t)(l - £L(t, T t )) . 

i=i 

For the sake of brevity, we suppress V YliLi n i * n ^ e above equation, which is irrelevant 
for hedging, and change the negative sign before the LIBOR futures to positive as follows 

N 

U(t) = CapfaU^ + V^THtyiLfaTi)) 

i=i 

N 

= Cap(t,U,T) + Vj2 n i(-t) f(t ' x) -l) ■ (20) 



3.1 Residual Variance 



Hedging a Cap denoted Cap(t, T) using LIBOR futures contracts can be accomplished by 
minimizing the residual variance of the hedged portfolio. It is the instantaneous change in 
the portfolio value that is stochastic. Therefore, the volatility of this change is computed 



to ascertain the efficacy of the hedge portfolio. 



The variance of the portfolio change, Var 



gift) 

dt 



equals 



Var 



dCap(t, T) 
Jt 



+ Var 



' N 

E 

i=l 



n 



dL{t,Tj t 
dt 



N 



riiVar 



i=i 



/dCap(t,T) dL(t,Ti] 
\ Jt ' Jt 



dCap(t,T) \ /dL(t,Ti] 



dt 



dt 



(21) 



The detailed calculation for determining the hedge parameters and portfolio variance 
is carried out in the Appendix. As in Baaquie, Srikant, and Warachka the following 
notation is introduced for simplicity 



M, 



rT+l rTi+£ 

xL{t,Tj) I dx dx'a{t,x)a(t,x')D{x,x';t) 
L(t,Ti)L(t,Tj) I dx J dx'a(t,x)a(t,x')D(x,x';t) 



(22) 



Equation (|22|) allows the residual variance in equation (|5()|) to be succinctly expressed as 
„T+£ r T+e N N N 



x 2 



/ dx dx'a(t,x)a(t,x')D(x,x';t) + 2 VAj^ + V V A i A j M ij (23) 

»=i i=i i=i 

which contains covariance terms. When at-the-money, the value of x below facilitates 
our empirical estimation of the model in Section 4 

r+oo -i / rT+e 2 

* = - VB(t ' T) L^A G -L dxmx) -T 



i -1/24, 



V 2 ^ 



^_ p -l/2^_ +iv(rf _ ; 



(24) 



v'w v 1 + ^ 

where d± = (in ~ ±q 2 /2j jq. The value of % for an at-the-money options yields d± = 
±g/2 which implies 

x(t,T)| at . the . money = VB{t,T)N{dJ) . (25) 

Observe that the residual variance depends on the correlation between forward rates 
described by the propagator. Ultimately, the effectiveness of the hedge portfolio is an 
empirical question since perfect hedging is not possible without shorting the original bond. 
This empirical question is addressed in Section 0] when the propagator is calibrated to 
market data. 



Hedge parameters rii that minimize the residual variance in equation ({23}) are 

N 

m = -Y^KjM^ 1 . (26) 

These parameters represent the optimal amounts of the futures contracts to include in 
the hedge portfolio. 

Equation (j2Ej) is proved by differentiating equation (|23*|l with respect to and sub- 
sequently solving for its value. The variance of the hedged portfolio in equation (j2"Tj) is 



proved by substituting the result of equation (|27)|) into equation (J23J) 

i-T+e pT+e N N 

V R = x 2 dx dx'a(t,x)a(t,x')D(x,x';t) -^^KiMr^Kj (27) 
•> T -> T i=\ j=i 

which declines monotonically as N increases. 

The residual variance in equation (|27|) enables the effectiveness of the hedge portfolio 
to be evaluated. Therefore, equation (|27jl is the basis for studying the impact of including 
different LIBOR futures contracts in the hedge portfolio. For N = 1, a single maturity 
Tj is evaluated, and the residual variance in equation (|2*7)l reduces to 

dx'a(t, x)a(t, x')D(x, x'\ t) 



X 



T+l r T+l 

1 dx 



T 



T 



^ ^ Jt +£ dx j^ +£ dx'cr(t, x)cr(t, x')D(x, x'\ t) \ ^ 
J^ +e dx J^ +£ dx'a(t, x)a(t, x')D(x, x'\ t) 



(28) 



The second term in equation represents the reduction in variance attributable to the 

hedge portfolio. To obtain the HJM limit, the propagator is constrained to equal one, 

reducing the residual variance Vr in equation (|2*H]l 

/ (J^ +e dx f^ 1+£ dx'cr(t,x)a(t,x')) 2 



X 



dxa(t, x) 



(29) 



it / \ f T * dxa(t,x)a(t,x') 

to zero. This HJM limit is consistent with our intuition that the residual variance is 



identical zero for any LIBOR maturity since all forward rates are perfectly correlated. 
This result is also shown empirically in SectionHJ However, results from hedging with two 
LIBOR futures contracts in HJM model are not presented since one degree of freedom 
cannot be hedged with two instruments. Indeed, in this circumstance, Mj" 1 is singular. 



3.2 Stochastic Hedging 



Residual variance enables us to control the effectiveness of the hedging procedure. How- 
ever, instead of only hedging downward movements in the Cap price, residual variance 
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operates on all forward rate fluctuations, including those that increase the portfolio's 
value. For this reason, we study stochastic hedging which is more practical since we can 
decide which forward rates to hedge against. 

Stochastic hedging of interest rate derivatives has been introduced by Baaquie [I], 
where the specific case of hedging Treasury Bonds is considered in detail. We focus 
on applying this technique to the hedging of a LIBOR Cap. Consider the hedging of 
a Cap against fluctuations in the forward rate f(t,x). A portfolio H(t) composed of a 
Cap(to,t*,T) and one LIBOR futures contract is studied. 

As in equation (J20J), we set N = 1 to obtain 



The portfolio is required to be independent of small changes in the forward rate. Thus, 
Delta hedging this portfolio requires 



In field theory, for each time t, there are infinitely many random variables driving forward 
rates, and one can never exactly Delta hedge by satisfying equation (J3*0*|) . The best 
alternative is to Delta hedge on average, and this scheme is referred to as stochastic Delta 
hedging as detailed in Baaquie P . To implement stochastic Delta hedging, one considers 
the conditional expectation value of the portfolio II (t), conditioned on the occurrence of 
some specific value of the forward rate fh = f(t,Xh), namely E[H(t)\f(t,Xh)]. Define the 
conditional probability of a Cap and a LIBOR futures by 
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U(t) = . 



(30) 



6f(t,x) 



Cap(t,U,T; f h ) 



E[Cap(t,U,T)\f h ] 



(31) 



L(t,T i; f h ) 



E[L(t,T x )\f h ] . 
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From Baaquie pQ and equation we have the conditional probability of a Cap 

given by 

/oo 
dG{(x-e G )+*(G\f h )} (32) 
-oo 

!k_ e - q h 2 e^G- q h J Dfe-it f^)^PiT +l ^f^)5{f{t,x h ) - f)e s 
= fDf6(f(t,x h )-f)e? ' 

while the conditional probability of a LIBOR futures is 



/oo 
dGe G ^{G\f ] t 1 T lj 
-oo 



/ - f£ +£ f(t, x)dx)S(f(t, x h ) - fY s 



* (G|/i( ' Tl) = WM-y ' (33) 



Stochastic Delta hedging is defined by approximating equation ("3T)j) as 

-§fEm)\h] = 0. (34) 

Hence, from equation ("3*""|) . stochastic Delta hedging yields 

dClp{t,U,T-,f h ) fiLjt^h) 

711 ~ W h 1 w h • ( } 

Thus, changes in the hedged portfolio Il(t) are, on average, sensitive to fluctuations in 
the forward rate f(t,Xh)- 

The conditional probability in equation ("3~2~|) and equation ([3*3)) along with the hedge 
parameter m is evaluated explicitly for the field theory description of forward rates in the 
Appendix which also contains the relevant notation. One should notice that nontrivial 
correlations appear in all the terms. The final result, from equation "55J1 . is given by 

C ■ C~ap(t,U,T; f h )-B-x-V ■ \xN'(d + )/Q + e - Go+a ^N(d_) - e~ G « +Q ^ N\d^) /Q 



ni = & 

As a comparison, the HJM limit is also analyzed in the Appendix. 
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Furthermore, one can Gamma hedge the same forward rate. To hedge against the 
d 2 U(t)/df 2 fluctuations, one needs to form a portfolio with two LIBOR futures contracts 
that minimizes the change in the value of £'[n(t)|// l ] by both Delta and Gamma hedging. 
These parameters are solved analytically, with empirical results presented in Section 4. 

Suppose a Cap needs to be hedged against the fluctuations of two forward rates, 
namely f(t, Xi) for i — 1, 2. The conditional probabilities for the Cap and LIBOR futures, 
with two forward rates fixed at f(t,Xi) = fa, are 

Cap(t,U,T-f u f 2 ) = E[Cap{t,U,T)\f u f 2 ] 
L(t,T i; A,/ 2 ) = ElLfrTJlhJz]. 

A portfolio of two LIBOR futures contracts with different maturities Tj 7^ T is defined as 

N 

IL(t) =Cap(t,t*,T) + J2m(t)Ht,Ti), (37) 

i=l 

where the hedging of this portfolio at instant time t is given by 

mt , fu A) = f* + ± "jff, + i ± ^f, + 1^ WA + o<a (38) 

with 5t = e = 1/360 year, while higher orders of e are negligible. Furthermore, the 
dynamics II = 5U/St equal 

HC, A.A) = f + t %h + g/^ + |^AA + O(e) ■ (39) 

Since (//) ~ i as in Baaquie [Tj, e/? ~ 0(1) ~ e/1/2, the second order terms are as 
important as the first order terms. Normal calculus retains the first order terms since e 
is infinitesimally small. However, e = 1 day in our context. 

The stochastic Delta hedging conditions are given by 

-^E[U(t)\f 1 J 2 ]=0 for j — 1,2 
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while stochastic Gamma hedging involves 

d 2 



Of 2 

with Cross Gamma hedging 

d 2 



2 E[n(t)|/ 1 ,/ 2 ]=0 for j = 1,2 



-£[n(t)|/i,/ 2 ]=o 



dfidh 

being unique to this paper. This Cross Gamma hedging only make sense in field theory 
models where movements in any specific forward rate can be hedged. 

One can solve the above system of N simultaneous equations to determine the N hedge 
parameters denoted n». The volatility of the hedged portfolio is reduced by increasing 
the number of forward interest rates being hedged. 

For this portfolio, we can analytically prove that Delta hedge parameters for the two 
forward rates differ by a prefactor 

^[n(0|/!,/ 2 ] = -^±E\a(t)\fuf2] = o (40) 

where A 2 and A 12 are defined in Appendix D. Therefore, Delta hedging against two for- 
ward rates can only determine the portfolio including one LIBOR futures. Furthermore, 
Gamma hedging two forward rates is the same except for a prefactor. 

Overall, for hedging against two forward rates we are left with three independent 
constraints from the above six constraints. In order to study the effect of each set of 
constraints separately, we form portfolios which include two LIBOR futures, and adopt 
hedging strategies that involve more than Delta hedging to fix the two hedge parameters. 
The first strategy implements one Delta hedge and one Gamma hedge on one forward 
rate. The two hedging parameters can also be fixed by one stochastic Delta hedge and 
an additional Cross Gamma hedge. 
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All of these hedge strategies are evaluated explicitly in the Appendix. Intuitively, we 
expect the portfolio to be hedged more effectively with the inclusion of the Cross Gamma 
parameter. Generally speaking, the field theory framework allows us to form portfolios 
that include more LIBOR futures and hedge against any number of forward rates. 

Until now, we obtained the parameter for each choice of the LIBOR futures and 
forward rates being hedged. Furthermore, we can minimize the following 

N 

Ew ( 41 ) 

i=l 

to find the minimum portfolio. This additional constraint finds the most effective futures 
contracts, where effectiveness is measured by requiring the smallest amount of contracts. 

In general, stochastic Delta hedging against N forward rates for large N is compli- 
cated, and closed-form solutions are difficult to obtain. 

4 Empirical Implementation 

This section illustrates the implementation of our field theory model and provides pre- 
liminary results for the impact of correlation on the hedge parameters. The correlation 
parameter for the propagator of LIBOR rates is estimated from historical data on LIBOR 
futures and at-the-money options. We calibrate the term structure of the volatility, cr(6), 
(see |H], [Oj) and the propagator with the parameters A and \i as in Baaquie and Bouchaud 
jHj. All the empirical results showed below are calculated from the derivation expressed 
in this paper. 
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4.1 Empirical Results on Residual Variance 

The reduction in variance achievable by hedging a Cap with LIBOR futures is the focus 
of this section. The portfolio 

N 

n(t) = Cap(t,U,T) + J2 n i( t )F( t , T i) 



i=l 



is considered with Var 



dU(t) 
dt 



being minimized. The residual variance for hedging a 1 
and 4 year Cap with a LIBOR futures is shown in Fig. Q along with its HJM counterpart. 
Observe that the residual variance drops to exactly zero when the same maturity LIBOR 
futures is used to hedge the Cap. 



3000 



2500- 



E 2000 



1500- 



1000 



500- 



- Hedging in HJM 
Hedging 1 year cap in FT model 
Hedging 4 year cap in FT model 



3 4 5 

Time to maturity/year 



Figure 1: Residual variance Var 



dU(t) 
dt 



for a one and four year Cap versus LIBOR 



futures maturity T x used to hedge portfolio U(t) = Cap(t,t*) + n 1 (t)J-'(t,Ti). 

By considering the changes of residual variance with respect of parameters A and fi, we 
find the neighboring points create no disparities, at least one cannot tell which offers the 
better hedge. An explanation of this is effect is that forward rates with similar maturities 
are strongly correlated. Furthermore, the HJM residual variance for both hedging a 1 
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Figure 2: Residual variance Var 



dU(t) 
dt 



for a four year Cap versus two LIBOR futures 



maturities Ti used to hedge U(t) = Cap(t,A) + Ylf=i n i{t)^{t, Ti). 

year and 4 year Cap are identical to the residual variance=0 axis. This is consistent with 
our analytical result in equation (J2HJ). 

The residual variance for hedging a 4 year Cap with two LIBOR futures is provided 
in Fig. |21 It is interesting to note that hedging with two instruments, even with similar 
maturities, entails a significant decrease in residual variance compared to hedging with 
one futures. This is illustrated in Fig. El where 9 = 6' represents hedging with one 
LIBOR futures. The residual variance in this situation is higher than the nearby points, 
and increases in a discontinuous manner. 



4.2 Empirical Results on Stochastic Hedging 

Stochastic hedging mitigates the risk of fluctuations in specified forward rates. The focus 
of this section is on the stochastic hedge parameters rij, with the best strategy chosen to 
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Figure 3: Hedge parameter ri\ for stochastic Delta hedging of Cap(t, 1,4) using LI- 
BOR futures maturity Ti and forward rate maturity Xh involving n(t) = Cap(t, 1,4) + 

ensure the LIBOR futures portfolio involves the smallest possible long and short positions 
since YliLi \ n i\ * s minimized. 

4.2.1 Hedging in Field Theory Models Compared to HJM 

The comparison is carried out in the simplest portfolio where one forward rate is hedged by 
one LIBOR futures, with a detailed empirical study in Subsection l4.2.2l As an illustration, 
Fig. El plots the hedge parameter n\ in our field theory model against the LIBOR futures 
maturity Ti, and the forward rate maturity x^ being hedged. One advantage of the field 
theory model is that, in principle, a hedge strategy against the movements of infinitely 
many correlated forward rates is available. To illustrate the contrast between our field 
theory model and a single-factor HJM model, we plot the identical hedge portfolio as 
above when D = 1, which has been shown to be the HJM limit of field theory models. 
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Figure 4: Hedge parameter n\ for stochastic hedging ofCap(t, 1,4) using LIBOR futures 
maturity T\ and forward rate maturity Xh in the HJM limit of D = 1 (forward rates 
perfectly correlated) involving n(t) = Cap(t, 1,4) + ni(t)J r (t,Ti) . 

From Fig. for the HJM limit, the hedge parameter n\ is invariant to the forward rate 
maturity Xh, which is expected since all forward rates f(t,Xh) are perfectly correlated in 
a single-factor HJM model. Therefore, it makes no difference which of the forward rates 
is being hedged. 

4.2.2 Hedging Against One Forward Rate with One LIBOR Futures 

We first study a portfolio with one LIBOR futures and one Cap to hedge against a single 
term structure movement. The portfolio is given by 

II(t) = Cap{t, t m ,T) + n^t)^, 7\) 

where the hedging is done by stochastic Delta hedging ■^■E\n.(t)\fh] = on forward rate 
f(t,x h ). 

Hedge parameters n\ for different LIBOR futures maturities 7\, and the forward rate 
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maturity Xh, are shown in Fig. El This figure describes the selection of the LIBOR futures 
in the minimum portfolio that requires the fewest number of long and short positions. 
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Figure 5: Hedge parameter n\ for stochastic hedging of Cap(t, 1,4) for forward maturity 
x h of forward rate f(t,Xh), with fixed LIBOR futures contract maturity T 1; involving 
II(t) = Cap(t, 1,4) + n 1 (t)F(t,T 1 ). 

Fig. El shows how the hedge parameters depend on Xh for a fixed T. Two limits 
T\ — 8 — j (3 months) and Ti = 165 are chosen. We find that Xh = 5 is always the 
most important forward rate to hedge against. Another graph describing the parameter 
dependence on 7\ is given in Fig. El with Xh = S. The minimum of hedge parameter 
ni at Xh — 1.5years reflects the maximum of o~(t,x) around the same future time. For 
greater generality, we also hedge Capit.U.T) for different and T values, and find that 
although the value of the parameter changes slightly, the shape of the parameter surface 
is almost identical. 
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Figure 6: Hedge parameter n± for stochastic hedging of Cap(t, 1,4) for LIBOR futures 
maturity T x when hedging against f(t, t + 5) with 5 = 3/12, involving IT(t) = Cap(t, 1, 4) + 

4.2.3 Hedging Against One Forward Rate with Two LIBOR Futures 

In Fig. [?l we investigate hedging one forward rate with two LIBOR futures by employing 
both Delta and Gamma hedging. The portfolio is given by 

2 

U{t) = Cap{t,U,T) + Y,^(t)F(t,Ti) 

i=l 

where stochastic Delta hedging -^E[H(t) \fi] = and stochastic Gamma hedging ^E[Il(t)\fi\ = 
are employed. 

From the previous case, we can hedge against f(t,5) in order to obtain a minimum 
portfolio involving the least amount of short and long positions. The diagonal reports that 
two LIBOR futures with the same maturity reduces to Delta hedging with one LIBOR 
futures. The data from which Fig. [7| is plotted illustrates that selling 38 contracts of 
L(t,t + 65) and buying 71 L(t,t + 5) contracts identifies the minimum portfolio. More 
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explicitly, the variables in the portfolio are given as 





T 2 


Xhl 


m 


n 2 


1.5 year 


0.25 year 


0.25 year 


-38 


71 




Figure 7: Summation of absolute hedge parameters \n\\ + \n 2 \ for two LIBOR futures, T\ 
and T 2 . The portfolio H(t) = Cap(t, 1, 4) + Y^=i ni{t)T{t, Tj) involves a stochastic hedge 
against one forward rate with both Delta and Gamma hedging. 

4.2.4 Hedging Against Two Forward Rates with Two LIBOR Futures 

In addition, we consider hedging fluctuations in two forward rates. Specifically, we study 
a portfolio comprised of two LIBOR futures and one Caplet n(t) = Cap(t,t*,T) + 
Si=i n i(*)-7 r (*j Ti) where the parameters rii are fixed by Delta hedging -^E[Il(t)\fi, / 2 ] = 
and and Cross Gamma hedging 9 £ 2 g ^ E[H(i)\fi, f 2 ) = 0. 

The result is displayed in Fig. |H1 where we hedge against two short maturity forward 
rates, such as f(t, 5) and f(t, 28). Again the data from which Fig. IHJis plotted illustrates 
that buying 45 contracts of L(t,t + 156) and selling 25 L(t, t + 35) contracts forms the 
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minimum portfolio. More explicitly, the variables in the portfolio are given as 
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n 2 
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Figure 8: Summation of absolute hedge parameters \rii\ + \n 2 \ for two LIBOR futures 
maturities, T% and T 2 . The portfolio U(t) = Cap(t, 1,4) + Y^i=i n i{p)^ : {pi r ^i) involves 
stochastic hedging against two forward rates, with both Delta and Cross Gamma hedging. 
Fig. and Fig. |H1 result from summing the absolute values of the hedge parameters 

(as in equation ([41))) which depend on the maturities of the LIBOR futures Tj. The 
corresponding empirical results are consistent with our earlier discussion. 2 

5 Conclusion 

LIBOR-based Caps and Floors are important financial instruments for managing interest 

rate risk. However, the multiple payoffs underlying these contracts complicates their 
2 If we choose the hedge portfolio by minimizing Yli=i n ii we ^ n( ^ that the minimum portfolio requires 

1500 contracts (long the short maturity and short their long maturity counterparts). 
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pricing as the LIBOR term structure dynamics are not perfectly correlated. A field theory 
model which allows for imperfect correlation between every LIBOR maturity overcomes 
this difficulty while maintaining model parsimony. 

Furthermore, hedge parameters for the field theory model are provided for risk man- 
agement applications. Although the field theory model implies an incomplete market 
since hedging cannot be conducted with an infinite number of interest rate dependent 
securities in practice, the correlation structure between LIBOR rates is exploited to min- 
imize risk. An empirical illustration demonstrates the implementation of our model. 
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A Residual Variance 



First, consider the variance of a Cap in the field theory model. Define the Delta of the 

n dCap(t,T) 

^ aP ' d f(t,x)dx' aS X 

r+oo jn i / rT+i 2 

* s - VB{t - T) L -t^A g -L ixf(t ' x) —2 



xle 2 " 2 V JT 2 J(X-e- G ) + }. (42) 
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The result in equation ()18|) for the Cap price implies that 



dCap(t, T) 
dt 



- v ' B(i ' r) L ^Ww{ G ~L M 

x < e 29 



2 



df(t,x) 



dt 



dx 



X\ I dxa(t,x) + / dxa{t, x)A(t 1 x) 
>t Jt 



(43) 

(44) 
(45) 



with E d ^§^ 

at 

variance equals 



dxa(t,x)j dt since E[A(t,x)] = 0. Therefore, the resulting 



dCap(t, T) 



-E 



dCap(t, T) 



X I dxa{t, x)A(t, x) . 



(46) 



With £(•) = - representing a delta function, squaring this expression and invoking the 
property that E[A{t,x)A{t,x')] = 5(0)D(x, x';t) = D<yX £ t — results in the instantaneous 
Cap price variance being 



Var 



dCap(t, T) 



-x 



dx 



dx'a(t, x)D(x, x'; t)a(t, x') 



(47) 



The quantity e signifies a small step forward in time. The underlying intuition is that 
we are converting a portfolio of futures contracts to one involving another function of 
LIBOR rates. Then, the instantaneous variance of a LIBOR portfolio is considered. For 
a LIBOR portfolio, fl(t) = V£J2? =1 \L{t,Ti), the following result holds, 



dU(t) 
dt 



E 



dtl{t) 
dt 



N rT t +£ 

J^AiL(t,Ti) / dxa(t,x)A(t 
i=i 



x 



(4* 



where L(t,Ti) = Ve^ f( t > x ) dx = T v T . +j ^ and 



Var 



dU(t) 
dt 



N N 



-VV A J A i L(t,T i )L(t,T i ) / dx dxa(t,x)D(x,x';t)a(t,x') . 

6 i=l j=l ' Jt * Jt 3 

(49) 
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The (residual) variance of the hedged portfolio 

11(f) = Cap(t, T) + J2 U 



N 



i=l 



is then computed in a straightforward manner. Equation implies the hedged portfo- 
lio's variance equals 

rT+i pT+i 

X 2 dx dx'a(t, x)a(t, x')D(x, x';t) 
Jt Jt 
N pT+e pTi+t 

+2x^2AiL(t,Ti) J dx J dx'a(t,x)a(t,x')D(x,x';t) (50) 

N N pT t +e pTj+e 

+ AjAjL(t, Tj)L(t, Tj) / dx I dx'a(t,x)cr(t,x')D(x,x';t) . 

1=1 j=l 

B Conditional Probability of Hedging One Forward 



Rate 



Using the results of the Gaussian models in Baaquie pQ, after a straightforward but 
tedious calculation, the following is derived from equations and 



*(G\f h ) 
$(G\f;t h ,T nl ) 



X 



1 



cxp 



cxp 



1 



2Q 2 
1 



(G - G ) 2 



2Q( 



2 (G — G1) 2 



(51) 
(52) 
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The notations are shown as follow 



{f Tn f 1 C f th c 

-J dxf(t ,x)-j a{t,x) + -E + —(f(t ,x h ) + J dta(t,x h ) - f- — ) 

d + = (lnx + G )/Q ; tZ_ = (In a; + G Q - Q 2 )/Q 

G = dxf(t ,x)-F--{f(t ,x h )-C-f + dta(t,x h )) + ^ 

rTni+i r Bi rt h 

d = dxf(t ,x)+ / a(t,x) — (f(t ,x h )- / dta(t,x h ) - f) 

JT nl JMz A Jt 

rth 

A = dta(t,x h ) 2 D(t,x h ,x h ;T FR ) 

Jto 



B = a(t,x h )D(t,x h ,x;T FR )a(t,x) 
J Mi 

Bx = a(t,x h )D(t,x h ,x;T FR )a(t,x) 

J Mi 

C = a(t,x h )D(t,x h ,x;T FR )a(t,x) 

J Mi 

D — a(t, x)D(t, x, x ; T FR )a(t, x ) 



q 2 = / a(t,x)D(t,x,x ;T FR )a(t,x 
JQ2+Q4 

E — I a(t,x)D(t,x,x ;T FR )a(t,x ) 



rth rTn rTn+i ^ ^ 

F — / dt dx dx a(t,x)D(t,x,x ;T FR )a(t,x ) . 

Jto Jt h JT n 

The domain of integration is given in Figs. E]andE3 It can be seen that the uncondi- 
tional probability distribution for the Cap and LIBOR futures yields volatilities q 2 and 
D respectively. Hence the conditional expectation reduces the volatility of Cap by 

B 2 

and by -4- for the LIBOR futures. This result is expected since the constraint imposed 
by the requirement of a conditional probability reduces the allowed fluctuations of the 
instruments. 
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It could be the case that there is a special maturity time Xh which causes the largest 
reduction in conditional variance. The answer is found by minimizing the conditional 



variance 

Cap(t h ,U,T n ;f h ) = X V(xN(d + ) - e - Go+ ^N(d.)) (53) 
L(t h ,T nl ;f h ) = e Gl+ ^. (54) 

Recall the hedging parameter is given by equation (J35J) . Using equation (jo^jl and setting 
to = t, t h = t + e, we get an (instantaneous) stochastic Delta hedge parameter rji (t) equal 
to 

C ■ Cap(t,U,T n ; f h ) - B ■ x • V ■ \xN'(d+)/Q + e~ Go+ ^N(d^) - e - G( >+^ iV'(cL)/Q 

2 m 



C HJM Limit of Hedging Function 

The HJM-limit of the hedging functions is analyzed for the specific exponential function 
considered by Jarrow and Turnbull JTT] 

a hjm (t,x) = a e^ x - t) , (56) 
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which sets the propagator D(t,x,x ;Tfr) equal to one. It can be shown that 

2 

2/3 

B = ^ e -frH( e -Pn_ e -PT n ^( e 2flt h _ e 2fit ) 

Ip 
ok 



(J — ^0 c -0x h f c -/3t h _ e -pT n \/ e 2pt h _ e Wto\ 

2/3 

2 

]j = a o r e -/3T n i+e c -(3T nl y^ c 2l3t h c 2/3t ^ 

2/3 3 

2 

£ _ _^<L( e -/3^ _ e -/3t h \2/ e 2pt h _ e 2f3t ^ 

2/3 3 

2 

F = _^2_( e -' 3 ' r «+ £ _ e -# r »)( e -0 r » - e"' 3 *' 1 )(e 2/3i ' 1 - e 2/3 *°) . 
2(3 3 

The exponential volatility function given in equation (|56|) has the remarkable property, 
similar to the case found for the hedging of Treasury Bonds in Baaquie jl] , that 

<fi(hjm) = D hjm -S^ = 0. (57) 

Hence, the conditional probability for the LIBOR futures is deterministic. Indeed, once 
the forward rate f\ is fixed, the following identity is valid 



Lhj m (th,T n i] f h ) = L(t h ,T nl ) . (58) 

In other words, for the volatility function in equation (|56jl . the LIBOR futures for the 
HJM model is exactly determined by one of the forward rates. 

However, the conditional probability for the Cap is not deterministic since the volatil- 
ity from th to £*, before the Cap's expiration, is not compensated for by fixing the forward 
rate. 
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D Conditional Probability of Hedging Two Forward 



Rates 



When hedging against two forward rates, equations (}3*2j) and (}3*3*j) imply we have the 
conditional probability of a Cap given by 



*(G|/i,/ 2 



TOO dp 

J — CO 27T 



jDf 



It h n Wh>*)jpI% M dxf{t h ,x) U 2 



and the conditional probability of LIBOR being 



/ Df6(G - f% +£ f(t h , x)dx) UU S(f(t h , Xi ) - f i)e ? . 
nG|/i,/ 2 ,T ni ) = FTTTTTl — ;T77, , TT^ J = 1,2 (60) 



which yields 

*(G\flj2) 

<$>{G\h,h,T nj ) 
under the following notation 



X 



1 



exp 



1 



2ttQ 2 



: exp 



2Q 2 
1 



(G-G ) 2 



2Qf j; 



J = 1,2 



(61) 
(62) 



V" = {l + £k)V 



l + £k 

X = exp j- / dxf(t ,x)- [ a{t,x) + ]-E + £^(i?i 2 - ^p)} 
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d+ = (In a; + G )/Q ; d_ = {lux + G - Q 2 )/Q 

fT n +£ r> 2 

G = / dxf(t ,x) -F- J1(R 12 -C 12 ) + \ 

JT n A12 1 



Q 



r>2 

2 _ „2 -°12 

il2 

(' 1 ii j t 



Gj = / dxf(t ,x)+ / a(t,x) j^-Rn J = 1 , 2 

JT n] JMj A 12 

Q) = Dj -?k j = 1,2 

J ^12 

Ri = f(t ,Xi)+ / dta(t,Xi) - fi £ = 1,2 



to 



i?12 — i?i 

^2 



Ai = dta(t,Xi) 2 D(t,Xi,Xi]T FR ) £ = 1,2 
A 12 = / dta(t,x 1 )D(t,x 1 ,x 2 ;T FR )a(t,x 2 ) 



A 12 = A, 



to 

An 



A 



2 

Bi = / a(t,Xi)D(t,Xi,x;T FR )a(t,x) £ = 1,2 
-B12 = -Bi 

^2 

i% = / a(t,Xi)D(t,Xi,x;T FR )a(t,x) £ = 1,2; j = 1, 

-Si2j = Bij T~^i i = 1 ) 2 ---;5 

A 2 

Ci = a(t,Xi)D(t,Xi,x;T FR )a(t,x) £ = 1,2 

J Mi 

C n A\ 2 „ 
12 — L'l "T - 1^2 

^2 



Dj = / a(t,x)D(t,x,x ;T FR )a(t,x) j = 1,2 
q 2 = a(t,x)D(t,x,x';T FR )a(t,x) 



E = a(t, x)D(t, x, x ;T FR )a(t, x ) 

F — / (it / do; / dx a(t, x)D(t, x, x ;T FR )a(t, x 

J to J t h JTn 

The domain of integration is given in Figs El and E3 
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Figure 9: Domain of integration M 1; M 2 and integration cube Q\, Q 2 , Qa where the x 
axis has the same limit as its corresponding x axis. 

Furthermore, an iV-fold constraint on the instruments would further reduce the vari- 
ance of the instruments 

Cap(t h ,U,T n ; hJ 2 ) = X V(xN(d + )-e- Go+ ^N(d.)) (64) 
L(t h ,T nj ;h,f 2 ) = e^ + ^. (65) 
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